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We calculate the amplitude Ti for forward doubly- virtual Compton scattering in heavy-baryon 
chiral perturbation theory, to fourth order in the chiral expansion and with the leading contribution 
of the 7NA form factor. This provides a model-independent expression for the amplitude in the 
low- momentum region, which is the dominant one for its contribution to the Lamb shift. It allows 
us to significantly reduce the theoretical uncertainty in the proton polarisability contributions to the 
Lamb shift in muonic hydrogen. We also stress the importance of consistency between the definitions 
of the Born and structure parts of the amplitude. Our result leaves no room for any effect large 
enough to explain the discrepancy between proton charge radii as determined from muonic and 
normal hydrogen. 

PACS numbers: 31.30.jr, 12.29.Fe, 13.60. Fz, 14.20.Dh 



I. INTRODUCTION 

The recent determination of the proton charge radius 
from the Lamb shift in muonic hydrogen [T] gives a value 
that differs by about 5 standard deviations from the CO- 
DATA value [2] and from the results of recent electron 
scattering experiments [3]. This has generated a flurry 
of activity attempting to understand the origin of the 
discrepancy and whether it could be a signal of physics 
beyond the Standard Model. Recent reviews of the situ- 
ation can be found in Refs. 0H7J-) 

However before concluding that new physics is re- 
quired, it is essential to examine carefully any possible 
conventional explanations. As discussed in the reviews, 
many of the contributions to the Lamb shift are theo- 
retically well constrained. One place where some theo- 
retical uncertainty remains is the contribution of proton 
structure to two-photon exchange, specifically through 
the polarisability of the proton. 

The energy shift of an S*-wave hydrogenic state due 
to two-photon exchange can be expressed in terms of 
the spin-averaged amplitude for forward doubly-virtual 
Compton scattering (V 2 CS) [SJ. This comprises two ten- 
sor structures: 

= (-g^ + ) Ti(w, Q 2 ) 

+ ^(/-^)(*-^*>(».<«. 

(1) 

where p and q are the four-momenta of the proton and 
photon, respectively, M is the nucleon mass, Q 2 = —q 2 , 
and ui = p-q/M. Dispersion relations [HI HO] can be used 
to estimate the amplitudes Ti,2(w,Q 2 ) from the corre- 
sponding structure functions measured in inelastic elec- 
tron scattering. These parts of the contribution of proton 



structure are well determined from the available data, 
with one important exception: the dispersion relation 
for T\ does not converge and so it requires a subtrac- 
tion. This is normally done at uj — 0, introducing a de- 
pendence on the unmeasured amplitude 21(0, Q 2 ). The 
slope of this term at Q 2 = is given by a low-energy 
theorem (LET) in terms of the magnetic polarisability 
of the proton, j3 [TTlfLl] . Otherwise its form is unknown 
and nucleon elastic form-factors are often used to model 
its Q 2 dependence [lUJ [T3] . The same amplitude is also 
needed in calculations of the electromagnetic self-energies 
of the nucleons [15l [16] and it has recently discussed in 
this context [IT] . 

This approach to determining the amplitude for for- 
ward V 2 CS has been questioned in several recent papers. 
Miller et al. [IB] have suggested that off-shell form-factors 
of the proton could generate new large polarisability con- 
tributions to V 2 CS. They also questioned the validity of 
the LET for the slope. Carlson and Vanderhaeghen [T5] 
estimated the Q 2 dependence from a two-pion exchange 
model motivated by chiral perturbation theory. They 
also split the amplitude Ti(o;,Q 2 ) into nucleon "pole" 
and "non-pole" terms, in contrast to Pachucki [TO], who 
separated out the full nucleon Born contribution obtained 
from the Dirac equation with on-shell form factors. Fi- 
nally, Hill and Paz [20] used an effective field theory, non- 
relativistic QED, to calculate ?i(0, Q 2 ), implicitly veri- 
fying the LET. They also questioned the use of on-shell 
form factors to get the Q 2 dependence of the non-pole 
Born terms in Ti(0, Q 2 ). They did not attempt to esti- 
mate the polarisability contribution but suggested that 
the uncertainty in two-photon exchange could have been 
significantly underestimated. More recently, Carlson and 
Vanderhaeghen |21) have examined the constraints on 
2i,2(w, Q 2 ) from real Compton scattering and shown that 
these are incompatible with the model of Miller et al. [IB] • 

To obtain a model-independent result for the form 
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of Ti(0,Q 2 ) at small Q 2 , we calculate Ti(0,Q 2 ) within 
the framework of chiral perturbation theory (ChPT) to 
fourth order. The polarisability contribution to the Lamb 
shift has previously been calculated using heavy-baryon 
chiral perturbation theory (HBChPT) at third order and 
found to be small [22] . However it is known that fourth- 
order contributions to the polarisabilities of the proton 
are larger than naively expected |23j . In addition, the 
A resonance is known to make an important contribu- 
tion to the magnetic polarisability [24l426j . and its po- 
tential significance for the Lamb shift has been stressed 
by Pineda [27j. We also note that ChPT results for 
spin-dependent V 2 CS show rapid dependence on Q 2 at 



low momenta [28 30 . If the spin- independent amplitude 
Ti(0,Q 2 ) were to show similar behaviour, this could be 
important for the Lamb shift. These results imply that it 
is important to check whether higher-order effects could 
alter previous conclusions. We do so here using the same 
ChPT approach as recently applied to the analysis of 
real Compton scattering [31] , treating nucleons to fourth- 
order in HBChPT and including effects of the A up to 
fifth order in the "5-counting" introduced by Pascalutsa 
and Phillips [55] . This gives us an estimate of the "mass" 
in the form factor for 7i(0, Q 2 ) which is determined con- 
sistently with the value for j3 obtained in Ref. [31] . 

In view of the questions that have been raised about 
the LET for 7i(0, Q 2 ), we re-examine its derivation to 
clarify some of the issues related to it and to the sepa- 
ration of the amplitude into Born and "structure" parts. 
We also comment on how the LET is embodied in non- 
relativistic effective field theories such as HBChPT. 

Finally, given the low-energy nature of the theory, it 
cannot give reliable estimates of the behaviour for larger 
momenta but we give arguments in favour of a l/Q 2 
fall-off at large Q 2 , as also obtained from the operator- 
product expansion |32j . We use this to estimate the con- 
tribution to the Lamb shift and to show that it is too 
small to explain the deduced discrepancy. 



II. LOW-ENERGY LIMIT 

Before calculating the higher-order contributions of 
proton structure to T± , we first re-examine the low-energy 
constraints on the proton-structure parts of the ampli- 
tude. In the low-momentum limit, the amplitudes for real 
and virtual Compton scattering can be related through 
the limited number of tensor structures that respect co- 
variance and gauge invariance, and are even under cross- 
ing (q O —q', v O [i). The invariant amplitudes mul- 
tiplying them should also be free from unphysical kine- 
matic singularities [33J 131] (see also: Refs. [H][T3]). Un- 
der these constraints, the spin-independent amplitude 
can be expanded in terms of five basis tensors: 



where p — (p+p')/2 is the average of the initial and final 
proton momenta, and we have taken the initial and final 
photons to have the same virtualitiesj^] These tensors 
may be chosen to be: 

C = -q' -qg^ +e<i ,v , 

C = q' ■ qp^ -q-P W + ^v v ) + (« ■ p) 2 iT, 
*T = «V + q'-q g"Y - <?W<? + «"V"), 

= q ■ p[2 q ■ p q 2 g^ -q-p (q»q u + q'Wl 

-q 2 (pW + q^f) + q' • q{?<? + q'^Pl] , 
t£" = q^p v + (q ■ p)Y"«" -q-P q 2 {P^q" + q'^P")- 

(3) 

Note that only four of the tensors are linearly indepen- 
dent, but there is no way to eliminate any one of them 
without introducing kinematic singularities in the ampli- 
tudes T l [m [T3J [3J] . 

In calculating the two-photon-exchange interaction be- 
tween a lepton and a proton, it is convenient to split the 
Compton amplitude into a Born piece, expressed entirely 
in terms of couplings of a single photon to an on-shell pro- 
ton, and a "structure" piece involving, for example, the 
polarisabilities of the proton. Conventionally, the Born 
amplitude is calculated from a Dirac equation for the 
proton with Dirac and Pauli form factors, Fp>(Q 2 ) and 
Fp{Q 2 ) respectively (often denoted F\ and F 2 although 
this invites confusion with inelastic structure functions). 
The Born contributions to the spin-independent Comp- 
ton amplitud^] obtained in this way have the forms 

rpB _ 2 (P 2 q' ■q-(q-p) 2 )Fp(2F D +F P ) 
a 6 M 3 (s — M 2 ){u — M 2 ) 

b 2 q' ■ q Fp (2F D + F P ) - AM 2 F 2 D - 2q 2 F D F P 



If = e 



T r B = 0, 



M 3 (s- M 2 ){u- M 2 ) 



rpB _ 2 F P (F D + Fp) 



M 3 (s- M 2 ){u- M 2 )' 



F B = 0. 



(4) 



r^HV 



7»(P 2 ,<7 2 ,9 ■ 



P,q'-q)tt 



(2) 



As required for the amplitudes corresponding to the ten- 
sors of Eq. ^3j, these have poles only at the on-shell 
points, s = (p + q) 2 = M 2 and u = (p ~ q') 2 = M 2 . 

The standard treatment of the Lamb shift [10] sub- 
tracts this Born amplitude to leave a residual, structure 
amplitude that is free of nucleon poles and can be esti- 
mated from inclusive-scattering structure functions with 



1 If the virtualities are different, there is also a sixth basis tensor. 

2 This has been defined by taking a covariant spin average, T^ v = 
Tr[(/ + M)T 3 M "(^ + M)]/(8M 2 ), similar to that used in Ref. [8]. 
Up to an overall factor, it corresponds to the spin-independent 
amplitude in the Breit frame, where there is a unique spin axis 
defined by the (parallel) initial and final momenta of the proton. 
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the aid of dispersion relations. This approach has been 
dismissed by Hill and Paz [20] as "sticking in form fac- 
tors" . However it provides a well-defined expression for 
the Q 2 dependence of the Born amplitude and, provided 
that the non-Born piece is defined consistently with it, 
there is no problem with this choice. Below, we comment 
further on the need for consistency between the choice of 
Born amplitude and the form of the non-Born piece. 

First, we look at the LETs satisfied by the non-Born 
amplitude defined in this way. This amplitude starts at 
second order in the photon momenta, with terms propor- 
tional to t a and tf,, the only tensors of order q 2 : 



47T 

M 2 



(a + f3)C + 0(q 4 ), (5) 



where the proton states are normalised so that the Thom- 
son term is (e 2 /M)e' ■ e. The contributions of proton 
structure at this order are contained in the two constants 
multiplying these tensors. They have been expressed in 
terms of the electric and magnetic polarisabilities of the 
proton, a and /3, which can be determined from real 
Compton scattering at low energies. (For details, see 
the review: Ref. [3T].) 

In the case of forward V 2 CS, the standard amplitudes 
m Eq. ffl can be related to those multiplying the five 
tensors of Eq. ^ by 

Tl = q 2 T a - (q ■ p) 2 T b - q 4 T c - 2q 2 (q ■ p) 2 T dl 
T 2 =M 2 (q 2 T b + q 4 T e ). (6) 

The residual amplitudes T, ; =Ti~Tf satisfy the LETs 

Ti(w,Q 2 ) = An Q 2 f3 + 4tt co 2 (a + P) + 0{q 4 ), (7) 
T 2 (uj,Q 2 ) =4TrQ 2 (a + f3) + 0(q 4 ). (8) 

The LET for Ti(0, Q 2 ) is important in the standard ap- 
proach to the Lamb shift [TU] since the dispersion relation 
for Ti requires a subtraction and the LET constrains the 
low-Q 2 limit of the subtraction term. 

Any introduction of off-shell form factors along the 
lines suggested by Miller et al. [18] should respect these 
LETs. The fact that their results for T x (0,(5 2 ) do not 
is more than just a numerical inconsistency with real 
Compton scattering [21] , it represents a violation of gen- 
eral principles. 

We now return to the question of consistency between 
Born and non-Born amplitudes. The factors of the Breit- 
frame photon energy q ■ p in the basis tensors of Eq. |3| 
mean that the Born contributions to the full Compton 
amplitude T^ v contain pieces without poles correspond- 
ing to on-shcll intermediate proton states. This can be 
seen clearly in the case of forward V 2 CS, where the Born 
terms have the forms 



2 ) = 



M 



Q 4 (F D (Q 2 ) + F P (Q 2 )Y 
Q4 _ 4M 2 w 2 



— F D (Q 2 ) 2 



T 2 B (cu,Q 2 ) 



4e 2 MQ 2 
Qi _ 4M 2 w 2 

x [F D {Q 2 ) 2 + {Q 2 /4M 2 )F P (Q 2 ) 2 ] 



(9) 



(10) 



and the second piece of T B contains no poles at u> = 
±Q 2 /2M. 

Various authors (for recent examples, see Refs. [TT1IT51 
|2"0] ) have suggested that, rather than splitting the Comp- 
ton amplitude into "Born" and "structure" pieces, it is 
better separate it into "pole" and "non-pole" pieces. In 
that approach, the second term of Eq. ([9| would count 
as part of the non-pole amplitude, along with the polar- 
isabilities and other structure contributions. However, as 
shown by Walker- Loud et al. |17J . this separation is not 
unique since the form of the pole piece depends on the 
choice of basis tensors. In the basis of Eq. ^ above, 
the Fd(Q 2 ) 2 term arises from a pole term in the am- 
plitude T b and, furthermore, it is required to keep the 
full Born amplitude free from kinematic singularities at 
Q 2 = 0. All of this implies that the separation into Born 
plus structure is the more natural one. 

The role of this term can be seen more clearly by ex- 
panding it in powers of Q 2 to get 



F D (Q 2 



= 1 - 



2M 2 



2 +O(Q 4 ), (11) 



where r c is the charge radius of the proton (defined in 
terms of the slope of the Sachs electric form factor) and 
k = Fp(0) is its anomalous magnetic moment. This 
shows that the leading term in this piece of T B gives 
rise to a term proportional to q^q^/q 2 in the scattering 
amplitude. In the full amplitude this cancels against a 
similar singular term arising from T B . 

The next term in Fp>(Q 2 ) 2 , of order Q 2 , generates a 
contribution to the scattering amplitude that remains fi- 
nite for real photons, Q 2 —> 0. Shifting this contribution 
from the Born amplitude to the structure part would al- 
ter the coefficient of Q 2 in T\ and so would correspond to 
a definition of the magnetic polarisability /3 that differed 
from the conventional one. Carlson and Vanderhaeghen 
[T§] keep only the pole pieces of the Born amplitude but 
the form they take for Ti(0, Q 2 ) is inconsistent with this 
since it satisfies an LET with conventional definition of 
p. The fact that the order-Q 2 term in T^O, Q 2 ) contains 
this piece from the expansion of Fp>{Q 2 ) 2 leads Walker- 
Loud et al. (17) to use the same definition for the struc- 
ture part of the V 2 CS amplitude as we advocate, albeit 
without the explicit justification presented here. 

In the context of the HBChPT calculation described in 
the next section, consistency between the Born and struc- 
ture amplitudes as usually defined means that we must 
subtract the full Born amplitude of Eq. ([9| expanded to 
the appropriate order in Q 2 . Specifically, in addition to 
all pole terms, we need to subtract the terms in the ex- 
pansion of Fp>(Q 2 ) 2 . To do this we need the Dirac form 
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factor calculated to third order in the chiral expansion, 
which is given in Ref. |35j . 

The LETs of Eqs. [§) are built into the effec- 
tive field theory for a nonrelativistic particle through 
the constraints of relativity, often expressed in terms of 
"reparametrisation invariance" [36] . The Lagrangian of 
any such theory contain three types of structure with two 
derivatives of the photon fields [37, 38 j. Two of these con- 
tribute to the polarisabilities a + j3 and ft. In contrast, 
the third has a coefficient that is fixed in terms of lower- 
order constants that correspond to the charge radius and 
anomalous magnetic moment. It generates the order-Q 2 
term of the non-pole piece of Tf discussed above. This is 
illustrated by the amplitude calculated in nonrelativistic 
QED, Eq. (11) of Ref. [50], which contains both the non- 
pole Born piece arising from the second term of Eq. (11) 
and the magnetic polarisability. 



III. CHIRAL PERTURBATION THEORY 

The specific amplitude we consider here is 21(0, Q 2 ) for 
forward V 2 CS, which appears as the subtraction term in 
the dispersion relation for T\{uj, Q 2 ) [TO] and so forms 
one of the least-well determined contributions of proton 
structure to the Lamb shift [TM20] . Its leading term, 
of order Q 2 , is fixed by the LET discussed in the previ- 
ous section but contributions from its higher-order terms 
vanish in the real-photon limit and so are subject to no 
such constraints. If these terms were to grow rapidly with 
Q 2 , they could lead to unexpectedly large structure con- 
tributions to the Lamb shift. They have been calculated 
by Nevado and Pineda [Hj at third-order where, like the 
polarisabilities j35] |39j , they are predictions that do not 
require any two-photon low-energy constants. However 
there are large contributions at fourth-order to both spin- 
independent and spin-dependent real Compton scatter- 
ing [23] . [28H30] . In addition, the effects of A resonance are 
important for the magnetic polarisability and could also 
play a significant role in the Lamb shift [57] . To examine 
whether these effects could be similarly important for the 
Lamb shift, we have calculated the Q 4 piece of forward 
V 2 CS treating nucleons to fourth order in HBChPT. We 
also include the leading contribution of the 7NA form 
factor, a term of fifth order in the 5 counting of Ref. [21] • 

The calculation of Ti(0,Q 2 ) has a very close corre- 
spondence to that of forward real Compton scattering 
[23 |39l 00] . To take advantage of this, we choose a gauge 
in which the polarisation vectors of the virtual photons 
are purely space-like (in the proton rest frame), so that 
the conditions p ■ e — 0, q ■ e — which were used in the 
real case still hold. With this choice, the diagrams that 
contribute are the same, and in many cases the resulting 
expressions can be obtained from the corresponding ones 
for real photons with the aid of some straightfoward sub- 
stitutions of kinematic variables. A representative subset 
of the fourth-order nucleon diagrams is shown in Fig. 1; 
the full set is given in Fig. 1 of Ref. [ID] . (The third-order 



diagrams can be found in Fig. 1 of Ref. |22j.) 

To determine the "structure" piece of the amplitude, 
Ti, we need to separate it from the Born terms. This sep- 
aration is somewhat subtle in HBChPT, and cannot just 
be done on a diagram by diagram basis [40]. Some dia- 
grams, such as Fig. 1(a), can immediately be identified as 
parts of the "magnetic" Born term (the first, "pole" term 
of Eq. ([9])). However there is also a one-nucleon-reducible 
diagram in which one photon couples to the nucleon via 
a magnetic-moment coupling and the other via a pion 
loop, Fig. 1(c). As discussed in Ref. [31], this contributes 
to both the magnetic Born and structure pieces of real 
Compton scattering. The structure part arises from the 
energy dependence of the pion loop which corresponds to 
an "off-shell form factor" of the sort proposed by Miller 
et al. [18] . but with a calculable form. This illustrates 
the fact that such off-shell dependences cannot be distin- 
guished from structure effects such as polarisabilities. 

In the case of space-like virtual photons considered 
here, particular care is needed with the one-nucleon re- 
ducible diagrams since the the leading propagator for the 
intermediate nucleon is 2M/Q 2 . This means that such 
diagrams contribute at one order lower than in the ex- 
pansion of real Compton scattering. For example, the 
diagram Fig. 1(c) just discussed contributes at third or- 
der, where it forms part of the magnetic Born piece of T\ . 
At fourth order there is a very similar graph but with a 
second-order insertion in the nucleon propagator within 
the loop. This contributes to both Born and "struc- 
ture" pieces of the amplitude. To this order, all other 
reducible diagrams with the enhanced 2M/Q 2 propaga- 
tor (diagrams like Fig. 1(a) but with higher-order ver- 
tices) contribute only to the magnetic Born term. 

The sum of all graphs (at third and fourth orders) has 
no constant term; the coefficient of Q 2 is the same as the 
loop contribution to j3 [25] plus a piece proportional to 
the third-order loop contribution to the derivative of the 
single-photon form-factor Fjj. The divergences in these 
loop integrals are renormalised by photon-nucleon sea- 
gulls, Fig. 1(e), arising from terms in the fourth-order chi- 
ral Lagrangian. These include a term with the structure 
F^F^ which is the same one that appears in the renor- 
malisation of the magnetic polarisability /?. This term is 



a combination of O 



(4) 



3,91,93,118 



in the minimal Lagrangian 



of Fettes et al. [38 . That Lagrangian also contains terms 
with the structure v^v" F\uF„ which contribute to a + f3 
and so are not relevant to forward V 2 CS. 

In addition, the fourth-order heavy-baryon effective 
Lagrangian contains terms with coefficients that are fixed 
by relativistic invariance in terms of coefficients from the 
lower-order Lagrangians. A combination of two of these 
is needed here, A40 and A52 of Ref. [38] . These have the 
structure [D M , F^^Dx + h.c, which does not correspond 
to either of the polarisabilities of the nucleon. But, like 
the similar cm term in the norelativistic QED Lagrangian 
of Manohar [37], the coefficient of this is given in terms 
of the second-order anomalous magnetic moment of the 
proton and the third-order coefficient that contributes to 



FIG. 1. Some diagrams contributing to Compton scattering in HBChPT: (a) a one-nucleon-reducible tree graph; (b-d) a 
representative selection of 7rN loop graphs that contribute at fourth order; (e) the fourth-order seagull contribution to /3; (f, 
g) A contributions at order S 3 and (h) at S 5 . In these graphs the circle, triangle and square denote vertices from the second-, 
third- and fourth-order Lagrangians, respectively. 



the charge radius. In the present context, this seagull 
diagram renormalises the loop contributions of order Q 2 
to F D (Q 2 ) 2 , the "non-pole" Born term. 

As well as fourth-order terms in the chiral expansion, 
we have also examined the effect of explicitly including 
the A resonance since this low-lying excitation can play 
an important role in the responses of a nucleon to ex- 
ternal fields, particularly in the magnetic polarisability 
[2~4T - [2l)] . Working within the 5 expansion introduced by 
Pascalutsa and Phillips [55], we find that the leading 
(order-<5 4 ) A-pole graph, Fig. 1(f), has the same form 
as the seagull contribution to /3 and so can simply be 
absorbed in the corresponding low-energy constant. In a 
similar way, the leading 7rA loops such as Fig. 1(g) can 
to a good approximation be absorbed into the the values 
of the low-energy constants appearing in the fourth-order 
tadpole diagrams, Fig. 1(d). The full expression differs 
from the pieces already included in the fourth-order re- 
sult by about 16%, which is well within the uncertainty 
introduced by the rather poorly-known second-order con- 
stants, and so we do not give it explicitly. 

The one A-pole contribution we do include arises from 
Fig. 1(h), where the 7rN loop generates a 7NA form fac- 
tor. This diagram is of fifth order in 5, which is higher 
than the other terms considered here, but it provides the 
leading contribution of this form factor to Ti(0, Q 2 ) at 
order Q A . We have therefore included it as an estimate 
of the likely influence of the A. 

After subtracting Tj 8 as discussed above, we arrive at 
the following result for T±; 

Ti(0,Q 2 ) 



= An/3Q 2 



e 2 gl 



12807r/2 TO 3 



'480Tr 2 f 2 ml 



4c, 



c 2 - 2c 3 



M, 



(4 + 5/i s 



e 9m 9 a Sttna 



Full expressions for the contributions to the amplitude, 
not expanded in powers of Q 2 , are given in the Appendix. 



In evaluating this numerically, we take g A = 1.27, f n = 
92.2 MeV, and fi s = 0.88. The constants Cj from the 
second-order Lagrangian have been determined from 7rN 
scattering and the resulting values, as quoted by Bernard 
[H], are c Y = -0.9±g;§, c 2 = 3.3 ± 0.2, c 3 = -4,7±\% 
all in GcV -1 . However a more tightly constrained value 
of C3 = —4.78 ± 0.2 GeV -1 has been obtained from pp 
scatttering [45] , and we use this here. 

For the couplings of the A, we take the values obtained 
by Pascalutsa et al. [12 H3] of g M = 2.9 for the magnetic 
7NA coupling constant and j IN4 = 1.425. These cou- 
plings provide a good description of real Compton scat- 
tering below the A peak [3T] but they do depend on the 
choice of Lagrangian. Nonetheless a very similar value for 
the product g A g TNA has also been found to give a good 
description of Compton scattering within a purely non- 
relativistic framework for the A [46] . despite ~ 30% dif- 
ferences in the separate couplings. A more serious source 
of uncertainty is the fact that these analyses are mainly 
sensitive to the couplings close to the A peak, whereas 
we need the values at zero energy. To account for this we 
allow for a 20% error on g A c^na 

If we express our result in the form of a form factor, 



T^Q^-AirfiQ 2 1 



91 



0(Q 4 ) 



its slope is given by 



P 



3 x 10" 4 fm 3 



(455 ± 32 MeV) 2 



(13) 



(14) 



The magnetic polarisability (3 has recently been de- 
termined from a fit to real Compton scattering within 
the same framework as used here, giving (3 = (3.1 ± 
0.5) x 10~ 4 fm 3 (including the statistical and Baldin 
Sum Rule errors only) |31j . Using this and the values 
for the Ci and A couplings above, with their errors, gives 
Mp = 460±50 MeV for the mass parameter in T i(0, Q 2 ). 
This has a size ~ 3m„-, which is "natural" for an effect 
with important contributions from pion loops. 

Without the fourth-order and A contributions to the 
slope, 455 MeV would be replaced by 588 MeV in 



() 



Eq. (14 1. In view of this change from 



third to fourth 
of higher-order 
the uncertainty 
However the A 
er source of un- 
the S expansion 
error associated 
take account of 



order, the residual error due to neglect 
terms is expected to be comparable to 
arising from the low-energy constants, 
contributions are likely to provide a largi 
certainty as the expansion parameter in 
is not particularly small. If we tri ple t he 
with of the order-Q 4 term in Eq. (12 1 to 
these uncertainties, then we get 

M B = 460 ± 100 ± 40 MeV, (15) 

where the first error is from the order-Q 4 term and the 
second is due to ft. 



IV. LAMB SHIFT 

The calculation outlined in the previous section leads 
to a model- independent result for the V 2 CS amplitude 
Ti(0, Q 2 ) up to order Q 4 . This constrains the subtraction 
term in the dispersion relation for Ti(w, Q 2 ) in a way 
that does not rely on the ad hoc use of the magnetic 
form factor of the proton as in Ref. [TO]. However it is 
not sufficient to estimate the polarisability contribution 
to the muonic Lamb shift, since that requires an integral 
over Q 2 . 

Our chiral calculation gives an expression for T±(0, Q 2 ) 
as a function of Q 2 which can be found in the Appendix. 
However, this cannot be used for large Q 2 since it con- 
tains terms (such as the fourth-order seagull) that do not 



vanish as Q 



Such terms were not present in the 



third-order calculation of Nevado and Pineda [22] , where 
the only photon-nucleon seagull is the one that gives the 
Thomson limit and so is part of the Born amplitude. 

Other information on the form of T\ (0, Q 2 ) comes from 
the large-Q 2 limit where the operator product expansion 
leads to 1/Q 2 for the dominant behaviour [32] (as noted 
in Ref. [30] ). This form can also be understood from the 
quark counting rules that apply to the partonic regime 
[47l 148] . In this regime there is no single-quark contri- 
bution to 7i(0, Q 2 ), as can be seen from the exact can- 
cellation of the two terms in Eq. (|9| at u> = for point- 
like fermions. Instead one-gluon exchange between two 
quarks gives the dominant contribution and the counting 
rules show that this falls off like 1/Q 2 . 

Here we estimate the contribution to the Lamb shift by 
matching our result smoothly onto the expected high-Q 2 
behaviour. To do this we write 



T!(0,Q 2 ) =^PQ 2 F p {Q 2 ), 
and take the dipole form 



Fp(Q 2 ) = 



1 



1 + Q 2 /2M B 



(16) 



(17) 



for the form factor. This matching could be carried out 
at Q 2 = 0, in which case M g is just the value Mg deter- 
mined from the slope of the ChPT form factor in the pre- 
vious section. However nonanalytic terms arising from 



pion loops can lead to non-dipole dependences of form 
factors on Q 2 in the region Q 2 < (3m v ) 2 (49j [50]. In 
Fig. 2 we show the full ChPT result for the form factor 
Fp(Q 2 ) compared to dipole forms with masses Mp = 460 
and 510 MeV. We see that the full form factor and the ex- 
trapolation from Q 2 = do deviate for Q 2 > 0.03 GeV 2 , 
although the differences lie well within the uncertainties 
of our calculation. In fact the ChPT result in the region 
from Q 2 ~ 0.1 to nearly 0.3 GeV 2 has a form that is very 
similar to a dipole with a mass of about 510 MeV. If we 
match onto the dipole form anywhere in this region we 
get a mass close to this value. In estimating the contri- 
bution to the Lamb shift we have therefore used Eq. ( 17) 
with 



Mg = 485 ± 100 ± 40 ± 25 MeV, 



(18) 



where the first two errors arise from the order-Q 4 term 
and /3, as discussed above, and the last is from the choice 
of matching point for the extrapolation. The shaded re- 
gion in Fig. 2 corresponds to our estimate of the the un- 
certainties associated with our calculation of the order- 
ed 4 term by showing dipole forms with masses in the 
range 385 to 585 MeV. 
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FIG. 2. The slope of the form factor for V 2 CS, plotted in the 
form [1 — Fp(Q 2 )]/Q 2 to emphasise variations at small Q 2 . 
The solid (red) curve shows the full result of ChPT calcula- 
tion. The dashed (blue) curve is the dipole form of Eq. (17) 
fitted to the slope at Q 2 = 0. The dot-dashed (black) curve 
is a dipole with a mass Mp = 510 MeV. The shaded region 
shows dipole forms with masses in range 385 to 585 MeV. 



The contribution of the subtraction term to the Lamb 



7 



shift is, from Eq. (34) of Ref. [TU]0 



a EM 0(O) 2 



47r m 



dQ 2 



ri(Q,Q 2 ) 
Q 2 




1 - 1 



where m is the lepton mass and 



m 2 



i 

8^ 



m M r , 



Mr, 



(19) 



(20) 



is the square of the 2S wave function at the origin. Al- 



though we evaluate this taking the dipole form of Eq. ( 17 ) 
for all Q 2 , we have also checked that its integral for 
Mp ~ 510 MeV agrees to better than 1% with that of 
the full form matched onto a dipole at Q 2 = 0.2 GcV 2 . 
Note that the leptonic factor strongly weights the low-Q 2 
region of this integral, with 90% coming from the region 
Q 2 < 0.3 GeV 2 . This means that matching onto the cor- 
rect form for larger values of Q 2 is not crucial, provided 
that the assumed form does not have a much longer tail 
than the dipole. 



Evaluating the integral using Eq. ( 17) for Ti(0, Q ) we 



get 



AE suh = 4.2 ± 1.0 fxeV, 



(21) 



for the contribution to the Lamb shift in muonic hydro- 
gen. The error here includes our estimate of the uncer- 
tainties associated with higher-order effects, as discussed 
in the previous section. However we note that largest 
single uncertainty comes from the value for f3, which ap- 
pears both as an overall factor in Ti(0, Q 2 ) and in the 
mass parameter Mp, as given by Eq. (14). 

Our result for A£ su b is more than twice that obtained 
by Pachucki [10]. Most of the difference is due to the 
modern value for /3, which is about twice the one used in 
that work. There is also a further ~ 10% increase from 
the use of the form factor with the correct asymptotic 
form. Our result is smaller that the one of Carlson and 
Vanderhaeghen [TU] who used a similar value for /3 as ours 
but whose model form factor leads to Ti(0, Q 2 ) growing 
like InQ 2 asymptotically. 

Finally, to estimate the total two-photon exchange 
contribution to the Lamb shift, we add in the "elas- 
tic" (Born) and dispersive contributions from Pachucki's 
treatment [THl EI]. We use the updated numerical re- 
sults from Carlson and Vanderhaeghen [19] and, for con- 
sistency with the definition of /3, we reinstate the "non- 
pole" Born piece. The latter gives about 4.8 /xeV for 



3 Note that our definition of Ti differs from that of Ref. [10] by a 
factor of e 2 M. 



the form factor of Ref. [52j and so is comparable in im- 
portance to the subtraction term. The main sources of 
uncertainty are /3 in our subtraction term above and the 
dependence of the elastic contribution on the choice of 
parametrisation of the proton form factors (which we 
take as ±1.6 /J.eV, based on the results in Ref. [19). The 
bottom line for the energy is 



AE 2l = -33 ± 2 /xeV, 



(22) 



which is not significantly different from the values of 
about —35 fj,eV obtained in Refs. [TOlQT] but which does 
not rely on a model for the form factor in T^O, Q 2 )- 



V. SUMMARY 

In this work, we have calculated the V 2 CS amplitude 
Ti(0,Q 2 ) at fourth order in HBChPT, and including 
contributions of the A resonance. Our result provides 
a model-independent constraint on the low-momentum 
form of the structure part of this amplitude. Since the 
contribution of this amplitude to the Lamb shift is domi- 
nated by the low-momentum region, this can significantly 
reduce the theoretical uncertainty in contribution of the 
polarisability of the proton to the Lamb shift, which has 
been the subject of some recent debate [T8Tf2"T] . 

When our amplitude is matched smoothly onto the 
high-Q 2 behaviour expected in the partonic regime [20] . 
we obtain a contribution to the Lamb shift that, perhaps 
unsurprisingly, is similar in magnitude to previous, more 
model-dependent determinations [101 US] ■ 

In reaching our result we have re-examined the LETs 
for V 2 CS and how they are embodied in effective field 
theories. The constraints of covariance, gauge invariance 
and crossing symmetry as well as lack of kinematic singu- 
larities mean that the only structure terms of order Q 2 
are given by the electric and magnetic polarisabilities. 
This is reflected by the appearance in fourth-order effec- 
tive Lagrangians of only two types of seagull term with 
adjustable coefficients, corresponding to the two polar- 
isabilities. These Lagrangians do contain other seagull 
terms at this order, but these all have cocfficcnts that 
are fixed by the requirements of covariance. In partic- 
ular, there is one that contributes to Ti(0, Q 2 ) at order 
Q 2 , but only through the "non-pole" piece of the Born 
amplitude. It is important keep this term in the Born 
amplitude, rather than subsuming it into the structure 
part [19) . as otherwise the magnetic polarisability in the 
LET Eq. |7| will not correspond to its conventional def- 
inition. 

Our results leave no room for any large additional po- 
larisability effect arising from off-shell form factors, as 
suggested by Miller et al. [IS]. The form factor appear- 
ing in Ti(0, Q 2 ) falls off with Q 2 and its mass parameter 
has a natural size for a pion-loop effect. We see no sign 
of any rapid growth at low Q 2 that could lead to a large 
contribution to the Lamb shift. 



There is thus no evidence that effects of proton polar- 
isability can explain the difference between the proton 
charge radius as recently determined from the muonic 
Lamb shift [1 and that obtained from other experiments 
[U [3J. The resolution of this puzzle must lie elsewhere, 
perhaps in re-analyses of the older experiments as sug- 
gested by Lorenz et al. [53 . 
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Appendix: Full amplitude 



Our result for the the V 2 CS amplitude as a full func- 
tion of Q 2 can be expressed in the form 



(3) 



T 



(4) 



(A) 



(A.l) 



where the third-order, fourth-order and A-pole contribu- 
tions are: 



T 



(3) _ e 2 g 2 A m n 



IQTrf 2 



. tan 



iQ 2 



47T 



mlQ 2 +0{Q q 



7687r 2 m w /2 



Ti 4) = 4TrS/3m z n Q 



2 ^2 



2"£%{ 12(C3 " 2Cl)(Q2 + 12) + 16C2( ° 2 + 3) + 3 lk ( 6 ~ (12Ms + 13) °' 



12 



: tanh 



Q 2 



Q 2 + 4 



12(c 3 -2 Cl )(Q 2 +4) + C2 (Q 2 + 4) 2 



(A.2) 



-3^(( Ms + l)Q 2 (Q 2 +4) 



(A.3) 



7^(A) _ 2c 2 2 2 

Jl ~ (M 2 -M 2 )(M & + M N ) m ^ < 



fj A Q n N A ' ^ 7T 

(Af A + M N ) 



.9m + .9a 



(Q 2 + 4) tan" , 
96tt/2^ 



2VQ 2 



(A.4) 



where Q 2 has been expressed in units of m 2 . Here 5/? 
denotes the sum of fourth-order loop and counterterm 

—(4) 

contributions to p. The subsequent terms in T 1 are of 
order Q A and higher, and hence contribute to the form 
factor that is needed for the Lamb shift. The value of 



(5/3 is independent of the choice of renormalisation scale 
and it is adjusted so that the complete term proportional 
to Q 2 in T\ gives the observed magnetic polarisability, 
as determined by a fit to Compton scattering data in 
Ref. [31]. 
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